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1. Introduction 

Single scale quantities in renormalizable quantum field theories, such as anomalous dimensions 
and massless Wilson coefficients, are most simply represented in terms of (finite) harmonic sums. 
This representation holds at least up to 3-loop order for massless Yang-Mills theories [1]. The 
corresponding Feynman-parameter integrals are such that nested harmonic sums [2, 3] appear in a 
natural way, working in Mellin space. They are written either as S- or as Z-sums, 

s, w= e "^y* 

N>h>ii...>i k >0 i, 

7 fAA- V nti[sign(m ; )]'' " 1 

^mi,...Mk\ ly ) ~ ,| m J 
N>ii>i 2 ...>ik>0 i[ 

Indeed a wide class of single scale quantities, including the anomalous dimensions and massless 
Wilson coefficients for unpolarized and polarized space- and time-like processes to 3-loop order, 
the Wilson coefficients for the Drell-Yan process and pseudoscalar and scalar Higgs boson produc- 
tion in hadron scattering in the heavy quark mass limit, as well as the soft- and virtual corrections to 
Bhabha scattering in the on-mass-shell scheme to 2-loop order [4] can be represented in this way. 
Single scale massive quantities at 2 loops, like the unpolarized and polarized heavy-flavor Wilson 
coefficients in the region Q 2 ^>m 2 - although the respective Feynman parameter-integrals exhibit 
a more involved structure - belong to this class too [5-8]. Finite harmonic sums obey algebraic, 
cf. [9], and structural relations [10]. The compact representations being obtained in Mellin space 
have to be mapped to momentum-fraction space to use the respective quantities in experimental 
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analyzes. The Mellin inversion requires the analytic continuation of the harmonic sums [11] w.r.t. 
the Mellin index N € C. 

Calculating massive operator matrix elements in Mellin space, which contribute to the heavy- 
flavor Wilson coefficients in deeply inelastic scattering, [6-8, 12], new types of infinite sums occur 
if compared to massless calculations. In the latter case, summation algorithms as summer [3], 
nestedsums [13] and Xsummer [14] may be used to calculate the respective sums, summer 
and Xsummer are based on FORM [15], while nestedsums is based on GiNaC. The new sums 
which emerged in [6-8, 12] can be calculated in different ways. In Ref. [6], we choose analytic 
methods together with summation. Another way consists in applying general summation algo- 
rithms in the field of computer algebra. The beginning was made by Gosper's telescoping algo- 
rithm [16] for hypergeometric terms 1 . Concerning practical applications, Zeilberger's extension of 
Gosper's algorithm to creative telescoping [17, 18] can be considered as the breakthrough in sym- 
bolic summation. The recent summation package Sigma [19] that can be executed in the computer 
algebra system Mathematica is going to open up completely new possibilities in symbolic sum- 
mation: Based on Karr's nZ-difference fields [20] and further refinements [21], the package con- 
tains summation algorithms [22] that allow to attack not only hypergeometric sums, like Gosper's 
and Zeilberger's algorithms, but sums involving also indefinite nested sums, like, e.g., S-sums or 
Z-sums dTT| ). 

In this note we report on using this package for the summation of a series of new sums and 
comment on different techniques, which may be used as well. 



2. The underlying summation principles 

In this section, we discuss methods of the summation package Sigma that are relevant to 
discover and prove the identities given below. Similar to hypergeometric summation [18], Sigma 
relies on the following summation techniques. 

Indefinite summation by telescoping: Given an indefinite sum S(a) = YX = q f {k) , find g(j) such 
that 

f(j)=8(j + l)-8U) (2- 1 ) 
holds within the summation range < j < a. Then by telescoping, we get 

S(a)=g(a + l)-g(0). 

Example. For the summand 

(2j + k+N + 2)j\k\(j + k+N)\ 



fU) 



(j + k+\)(j + N+\)(j + k+l)\(j + N+\)\(k + N+l 



j\k\(j + k + N)\(-Si (j) + Si (j + k) + S 1 (j + N )-S 1 (j + k + N)) 

{j + k+l)\(j + N+l)l(k + N+l)\ { ' 



S i gma finds 

U + k + l){j+N+ l)j\k\(j + k + N)\{s l (j) - Si (j + k)-Si (j + N) + Si (j + k + N)^j 



kN(j + k+l)\(j + N+l)\(k+N+l) 



l f(k) is hypergeometric in k iff f(k + l)/f(k) = g{k) for some fixed rational function g(k). 
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Therefore summing (2. 1 ) over j produces 



£/Q) 



(a + 1) \{k - 1) \{a + k + N + 1) ! (Si (a) - S l (a + k) - Si (a + N) + Si (a + k + N)) 



p J w/ N(a + k+l)\(a + N+\)\{k + N+l)\ 

Si(k)+Si(N)-Si(k + N) 



+- 



kN(k+N+l)N\ 



{2a + k + N + 2)a\k\{a + k + N)\ 

{a + k+\)(a + N +\)(a + k+\)\(a + N +\)\{k + N +l)V 

for a G N. 

More generally, given a finite set of indefinite nested sums and products (in particular, hyper- 
geometric terms and Z-sums/S-sums) and given f(j) as a rational expression in terms of those 
sums and products, Sigma decides algorithmically, whether there exists a rational expression g(j) 



in terms of those sums and products such that (2.1) holds; for more details see [19]. 



Usually, one fails to find a solution g(j) for (2.1). If f(j), depends on an extra discrete param 



eter, say Af, Zeilberger's creative telescoping paradigm can be applied. 

Deriving recurrences by creative telescoping: Given an integer d > and given a sum 

S(a,N):=j^f(N,j) (2.4) 

with an extra parameter N ', find constants co(N), . . . ,Cd(N), free of j, and g(N,j) such that the 
following summand recurrence 

co(N)f(N, ;) + ..., c d (N)f(N + d, j) = g(N , ; + 1 ) - g(N, j) (2.5) 

holds for the summation range < j < a. If one succeeds in this task, one gets by telescoping the 
recurrence relation 

c (N)S(a,N) + --- + c d (N)S(a,N + d)=g(N,a + \)-g{N,0). 

Summarizing, one arrives at a recurrence relation of the form 

c Q (N)S(a,N) + ...c d (N)S(a,N + d) =q(a,N) (2.6) 

for the sum S(a,N) of order d. Note that a can be specialized, e.g., to Af, i.e., one obtains a 
recurrence for the definite sum Lf=o/(^^)- In other instances, see below, one can perform the 
limit a — > oo which yields a recurrences for Y,T=of(N,k). 

Creative telescoping has been introduced for hypergeometric terms [17]. More generally, in 
Sigma the summand f(N,j) may consist of indefinite nested sums and products, in particular 
hypergeometric terms and S-sums/Z-sums. 

Example. For d = 1 and the summand 

JK ,JI jN(j+N+l)N\ K ' 
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S i gma computes the solution 



c (N) = -N{N+ l) 2 , Cl (N)(N+ l) 3 {N + 2), 

jS(l,j) + (-N-l)S(\,N)-jS(l,j + N)-2 (2-8) 



g(NJ 



(j + N+\)N\ 



Hence, summing ( |2.5| ) over 1 < j < b gives 



■NS(N,b) + (1 + N)(2 + N)S(N+ 1,*) 



+ 1) (S, (b) + S, (N) - Sl (b + N)) + b(b+l) (2 9) 



(N + l) 2 (b + N + 2)N\ (N + iy(b + N+ l)(b + N + 2)N 

for the sum S(N,b). 

Summarizing, if we succeed in finding a recurrence 2 of type 

c (N)S(N) + ...c d (N)S{N + d)=q(N), (2.10) 

then together with the initial values S(i), < i < d, we obtain an alternative representation of the 
corresponding input sum, F- =0 (resp. l!J =Q f{NJ) or ZJ =Q f(N,j)). 

For many applications, such a result is completely satisfactory. However, if one hunts for a 
closed form, one can continue as follows. 



Recurrence solving: Given a recurrence §L.\0 ),find linearly independent solutions of the homo- 



geneous version of the recurrence, say h\(N), . . . ,h r (N), and a particular solution of (16), say 
p(N). 

Namely, if we manage to compute sufficiently many solutions, we can find constants ki, . . . ,k r 
such that 

S(i) = p(i) +k\h\(i) H \-k r h r (i) , 

for all < i < d. As a consequence, 

S(N) = p(N) +k\h\(N) H \-k r h r (N), 

for all N>0. 



With Sigma we can handle the following situation: Given ( |2.6| ) where the cq(N), . . . ,Cd(N) 
and q(N) are given by indefinite nested sums and products, Sigma can compute all solutions S(N) 
in terms of indefinite nested sums and products (the so-called d'Alembertian solution; see [23] 
and [19, Section 7.2.3]). We emphasize that Sigma finds, as a special case, all solutions in terms 
of hypergeometric expressions and 5-sums/Z-sums. 



3. Summation through difference equations: single sums 

As outlined in Refs. [6,7], in course of the calculation of massive 2-loop integrals containing as 
single (general) scale the Mellin parameter N, different types of sums occur. These summands are 

2 For simplicity, we suppress extra parameters, like e.g., a. Moreover, we assume that Cj(N) ^ for all WgN. 
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typically products of harmonic sums with different arguments, weighted by summation parameters 
and contain hypergeometric terms, like binomial factors or Beta-function factors 

Here i is the summation-index, which runs from one to infinity. In general, these sums can be 
expressed in terms of nested harmonic sums, [2, 3] and ^-values. Note that sums containing Beta- 
functions with different arguments, e.g. B(i,i), B(N + usually do not lead to harmonic sums 
in the final result. Further quite often binomial sums of the type 

iHf(i), (3-2) 

i=o V 1 J 

with F(f) being a broken-rational function of i times a product of harmonic sums, emerge. This 
type of sums can as well be treated with the Sigma-package but are not considered here. 

Some of these sums can be performed by the existing packages [3, 13, 14]. However, there 
exists so far no automatic computer program to calculate sums which contain Beta-function factors 
of the type (3T) and single harmonic sums in the summand. 

3.1 The Sigma-approach 

As a first example we consider the sum 

Ti (*) := t r^fi Si + o ■ (3 - 3) 

We treat the upper bound of the sum as a finite integer, i.e., we consider the truncated version 

" B(N,i) ,. , 
71 (a, AT) :=E . + y +2 fr(0ft(# + 0, 

for a G N. Given this sum as input, we execute Sigma's creative telescoping algorithm and find 
a recurrence for T\(a,N)(= S(a,N)) of the form (2.6) with order d = 4. Finally, we perform the 
limit a — > oo and we end up at the recurrence 

- N(N + 1 ) (N + 2) 2 (4N 5 + 6SN 4 + 455N 3 + 1494N 2 + 2402N + 1510) Ti (N) 
- (N + 1 ) (N + 2) (N + 3) ( 1 6N 5 + 26(W 4 + 1 660N 3 + 5 1 SSN 2 + 79 12N + 4699) T^N+l) 
+ (N + 2)(N + 4) (2N + 5) (4N b + 14N 5 + 542N 4 + 1978N 3 + 3680N 2 + 3 103N + 767) 7\ (N + 2) 
+ (N+4)(N + 5) (16N 6 + 276N 5 + \92SN 4 + 6968N 3 + 137 16N 2 + 13929N + 5707) 7\ (N + 3) 
- (N + 4) (JV + 5 ) 2 (A 7 + 6) (4N 5 + 48N 4 + 223N 3 + 491 N 2 + 527N + 2 1 1 ) 7\ + 4) 

= P l (N)+P 2 (N)S l (N) 

where 

P x (AT) = (32N 18 + 1232N 17 + 21512N 16 + 223472N 15 + 1514464N 14 + 68061 14N 13 

+ 18666770N 12 + 15297623N 11 - 116877645Af 10 - 6414589 BN 9 - 1826931522^ 
- 350720529 IN 7 - 4825457477N 6 - 4839106893N 5 - 3535231014N 4 - 1860247616N 3 

_ 684064448N 2 - 160164480N - 1739520o) / (N 3 (N + 1 ) 3 (N + 2) 3 (N + 3) 2 (N + 4) (N + 5)) 
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and 

P 2 (N) = -4^(4N U + 15(W 13 +261(W 12 + 27717/V 11 + 199197/V 10 + 1017704/V 9 + 3786588/V 8 

+ 10355813/V 7 +20779613/V 6 + 30225025/V 5 + 31132328/V 4 + 21872237/V 3 + 9912442/V 2 

+ 2672360/V + 362400) / (N 2 (N + l) 2 (N + 2) 2 (N + 3) (N + 4) (TV + 5)) . 

In the next step, we apply Sigma's recurrence solver to the computed recurrence and find the four 
linearly independent solutions 

1 (-If 



hi{N) 



h 3 (N) 



N + 2' 
Si (AO 



h 2 (N) 



h 4 (N) 



N(N +l)(/V + 2)' 

(l + (A r + i) 5l(A r))(_i) 



7V + 2' v ; 7V(/V+l) 2 (/V + 2) 

of the homogeneous version of the recurrence plus the particular solution 

2{-l) N 



p(N) 



N(N+l)(N + 2) 

2S_ 2 (AQ + C2 
^ 3 N+l 



25_2,i (iV) -35_3 (iV) -25_2 (N)5i(iV) - &5i(iV) 



N + 2 



N + 2 



Si(N) 



2 + lN + lN 2 + 5N 3 +N 4 2 + 1N + 9N 2 + 4N 3 +N 4 

+ N 3 (N+l) 3 (N + 2) l( )+ N 4 (N+l) 3 (N + 2) 
of the recurrence itself. Finally, we look for constants c\ , . . . ,c 4 such that 

T\ (N) = c\ h\ (N) + c 2 h 2 (N) + c 3 h 3 (N) + +c 4 h 4 (N) +p(N). 

The calculation of the necessary initial values for TV = 0, 1 , 2, 3 does not pose a problem for S i gma 
and we conclude that c\ = c 2 = c 3 = c 4 = 0. Hence the final result then reads 

2{-l) N 



7i(A0 



N(N+l)(N + 2) 

2S_ 2 (AQ + C 2 
^ "' N+l 



2S_ 2) i (N) - 3S- 3 (N) - 2S^ 2 (N)S l (N) - fcSj (N) 



M N )~b s 2 (N)-i; 2 



N + 2 



Si(N) 



2 + 7N + 7N 2 + 5N 3 +N* „ 
+ — S x (N) +2 



N + 2 

2 + 1N + 9N 2 + 4N 3 + N 4 



(3.4) 



N 3 (N+l) 3 (N + 2) 1V 1 N 4 (N+l) 3 (N + 2) 

Using more refined algorithms of Sigma, see e.g., [24], even a first order difference equation 
can be obtained 

(N + 2)7i (AO - (N + 3)7i (N + 1) 



= 2 



(-1) 



N 



3N + 4 



N(N + 2)y (N+l)(N + 2) 

S 2 (N) - £ 2 N 6 + SN 5 + 3 IN 4 + 66N 3 + 88/V 2 + 64N + 16 



■ + ■ 



+2 



N+l N 3 (N+l) 2 (N + 2) 3 

N 5 + 5N 4 + 21N 3 + 3SN 2 + 2SN + 8 
N 4 (N+l) 2 (N + 2) 2 ' 



£ 2 + 2S_ 2 (iV)) " 2^3 - 2S- 3 (N) - 2C 2 5i (AT) - 45 lj _ 2 (7Y)J 

Si (AO 

(3.5) 
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However, in setting up Eq. (3.5), use had to be made of further sums of less complexity, which 
had to be calculated separately. As above, we can easily solve the recurrence and obtain again the 



result (3.4). 



Here and in the following we applied various algebraic relations between harmonic sums to 
obtain a simplification of our results, cf. [9]. 

3.2 Alternative approaches 

As a second example we consider the sum 



r 2 (A0:=L 



S 2 (i + N) 



(3.6) 



which does not contain a Beta-function. In a first attempt, we proceed as with the first example 
Ti(N). Namely, the naive application of Sigma yields a fifth order difference equation, which is 
clearly too complex for this sum. However, similar to the situation T\ (N), Sigma can reduce it to 
a third order relation which reads 

T 2 (N) (N+ l) 2 — T 2 (N + 1) (3N 2 + 1(W + 9) 
+T 2 {N + 2) (3N 2 + 14N + 17) - T 2 (N + 3) (N + 3) 2 



6/V 5 + 48/V 4 + 143/V 3 + 186/V 2 + 81/V-12 2N 2 + 1N + 1 



(# + 1)2 (at + 2)3 (AT + 3)2 (^+1)^ + 2)2(^ + 3) 

2N 6 - 24N 5 -U6N 4 - 2SSN 3 - 386N 2 - 264N - 72 

(N+\) 2 (N + 2y(N + 3) 2 ^ ' 



Si (AO 



(3.7) 



Solving this recurrence relation in terms of harmonic sums gives a closed form; see (3.14) below. 



Still (3.7) represents a rather involved way to solve the problem. A better way consists in first 
mapping the numerator S 2 (i+N) into a linear representation, which can be achieved using Euler's 
relation 



S 2 a (N)=2S aja (N)-S 2a (N), a>0. 

This is realized in summer by the basis-command for general-type harmonic sums, 

2SiA{i + N)-S 2 (i + N) 



(=i 



r 



As outlined in Ref. [3], sums of this type can be evaluated by considering the difference 

1 

; __ ] — — -L — 



D^ = r 2 (y)-r 2 a-i) = 2£^-| i , ( . +02 . 



The solution is then obtained by summing the difference ( |3.10| ) to 



(3.8) 



(3.9) 



(3.10) 



T 2 (N) = ^D 2 (j) + T 2 (0) 
7=1 



(3.11) 
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The sums in Eq. (3.10) are now calculable trivially or are of less complexity than the original sum. 
In the case considered here, only the first sum on the left hand side is not trivial. However, after 
partial fractioning, one can repeat the same procedure, resulting into another difference equation, 



which is now easily calculable. Thus using this technique, the solution of Eq. ( |3.6[ ) can be obtained 
by summing two first order difference equations or solving a second order one. The above proce- 
dure is well known and some of the summation-algorithms of summer, [3], are based on it. As a 
consequence, infinite sums with an arbitrary number of harmonic sums with the same argument can 
be performed using this package. Note that sums containing harmonic sums with different argu- 
ments, see e.g Eq. ( |3.6| , |A.2| ), can in principle be summed automatically using the same approach. 
However, this feature is not yet built into summer. 



A third way to obtain the sum ( p.q ) consists of using integral representations of harmonic 
sums, [2]. One may represent the sum in terms of the following integrals 



1 J+N 



ln(l-x) 



ln(l-jt) 



W=2£/ n dx- 2 , l v 

Li 2 (x) (N+l) 



2M 



l-x 





l-x 



U z (x) (#+!) + £ . (3.12) 



Here the Mellin-transform is defined as 

M 



[f(x)](N) := tdx^fix) 
Jo 



(3.13) 



Eq. (3.12) can then be easily calculated since the corresponding Mellin-transforms are well- 



known, cf. [2]. Either of these three methods above leads to 

72 (AO = -^2 + 45i (A% + Si (AO £2 - S 2 (N)^ 2 - 25i (N)S 2A (N) - S 2 , 2 (N) . (3.14) 



As a third example we would like to evaluate the sum 

i=l { 



(3.15) 



Note that Eq. ( |3.15[ ) is divergent. In order to treat the divergent pieces, we introduce the symbol Ci 

(3.16) 



a 1 

Oi := lim V" — . 

a — ^00 7 



The application of Sigma to the sum ( 3.15| ) yields a fourth order difference equation 



(N+\f{N + 2)T 2 (N)-(N + 2) (4N 2 + \5N + 15) T 2 (N + 1) 

+(2N + 5) (3A^ 2 + 15^V + 20) T 2 (N + 2)-(N + 3) (4N 2 + 25A^ + 40) T 2 (N + 3) 

+(N + 3)(N + 4) 2 T 2 (N + 4) 

6N 5 + 13N 4 + 329N 3 + 6S4N 2 + 645N + 215 6N 2 + l9N + 9 

+ 7TTT^~~ — Zz Sl[N) , (3.17) 



(N+l) 2 {N + 2) 2 (N + 3) 2 



(N+l)(N + 2)(N + 3) 
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which can be solved. As in the foregoing example the better way to calculate the sum is to first 
change Sj(i + N) into a linear basis representation, cf. [3], 

Tm = i 2S ' Ai+N) .- Hi+N) S,(0- (3-18) 

i=\ 1 

One may now calculate T3 (N) using telescoping for the difference 



(3.19) 



with 



N 

T 3 (N) = ^D 2 (j) + T 3 (0) . (3.20) 
7=1 



One finally obtains 

m) _ ^ + gg +H , W 6 + M?(w) 2 - ftW fc-2s,Wj t ,W 

+5?(Ar)5 2 (Ar)+5 1 (Ar)5 3 (Ar)-SM + 3M . (3.21) 
4. A double sum example 

Using the Mellin-Barnes integral representation [25], massive two-loop integrals which occur in 
the calculation of polarized and unpolarized massive operator matrix elements [6-8] result into 
double infinite series of the kind 



00 00 



S(A0 = II/(Ar,M- (4-D 

k=0j=0 

Here N € N denotes the Mellin variable. In the following, we consider an example dealt with 
in [6, 8] where 

f(N k j) = F{k+l) r(e)T(l - e fU+l-2e)rU+l + emk + j+ l+N) 
J[?i,k,J) r( k + 2+ N) [ } { ' T{j+l-e)T{j + 2 + N)T{k + j + 2) 

T(k+l) T{j+ i + 2e)T{j+l-e)T(k + j+l + e + N) 

^T{k + 2 + N) y ] K ^ ' r{j+l)rU + 2 + e + N)r{k + j + 2 + e) ' 

e denoting the dimensional regularization parameter of the momentum integrals performed in 
D = 4 + e space-time dimensions. f(N,k,j) obeys the expansion 

f(N,kJ)=MN,k,j) + eMN,k,j) + e 2 MN,k, + , (4.2) 

which is derived using 

w^r l+ % etz -'' (n - ,h <43) 
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with It = {1,...,1}, cf. e.g. [2,26]. Since the summation could not be accomplished with the 

k 

available packages summer [3], nestedsums [13], or Xsummer [14], this double sum was 
solved in [6, 8] using suitable integral representations, which allowed the summation. Here special 
partial differential operators had to be used to map the series expansion efficiently. This method 
clearly was specialized to the respective cases to be dealt with. 

Let us consider the summation of the leading term fo(N,k,j). The summation for the higher 
terms in e can be performed in an analogous way. 

We start with the summand fo(N,kJ) which is given in (2.2), and find the closed form given 



in (|2.3D. 

Then the limit 

is performed leaving us with an expression in finite harmonic sums which depend on the second 
summation parameter k. We repeat the foregoing procedure and form a finite sum 

b G N. 

First, we execute Sigma's creative telescoping algorithm and compute the recurrence rela- 
tion (2.9). Again, by limit considerations, it follows that the sum 



oo oo 



S'(N) := lim S(N,b) = £ £ fo(N,k,j) (4.6) 
^- k=ij-o 

satisfies the recurrence relation 

-NS'(N) + (1 +N)(2 + N)S\N+ 1) = ^+i^±l . 

Next, we use Sigma's recurrence solver and compute the general solution 

1 S l (N) 2 + S 2 (N) 

c+ ■ 



Af(Af+l)AH 2N(N+l)N\ 
for a constant c. Checking the initial value S'(l) shows that 

S l (N) 2 + S 2 (N) 
[ V 2N(N+l)l ' 

Hence, 
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Using again S i gma we find 



Y f (n n ■\ = V V+N + 2 

LM»,»,J) L u+mj+N+mN+1)l 



7=0 j=0 

2a + N + 2 



(a+l) 2 (a + N+l) 2 (N+l)\ 

S 2 {a) S 2 (N) S 2 (a + N) 

N(N+l)\ N(N+\)\ N(N+1)\' 



In the limit a — > oo this later expression simplifies to 

S 2 (N) 



N(N+1)\ 

Finally, we find 

Si(A0 2 + 3S 2 (A0 



Completely analogously, we compute, for instance, the linear term in the series expansion (4.2): 



ff fftf , , S l (N)*-3S 2 (N)S l (N)-ZS 3 (N) 

hg> = ^Ti)! • 

A. Appendix 

Further examples of sums calculable by S igma and the order of their respective recurrence relation 
are given below. Note that the inhomogeneous part of the recurrence relations are given in terms of 
single harmonic sums up to depth 2, ^-values and polynomials in N only. Below we list a series 
of examples for infinite sums, which emerged calculating the linear terms in e for the massive 
operator matrix elements [12]. The sums were calculated solving higher order difference equations 
with Sigma. 

A.l Sums resulting in fourth order difference equations 



i=l 1 1U 

+2s^ 2 (iv) - 45, (mm - siwhm + ^ 
= ^ 6 _»g_ SlCW)6+ > fecw) _ s?(w)1Cl+ ^ cw) 



i=l 



+Sj(N)S 2 (N) + S l (N)S 3 (N) . (A3) 
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A.2 Sums resulting in third order difference equations 



y v 



g(AM) 

i + 2) 



Mi) 



A , 25i,- 2 (jV)+5-3(jV) + C 2 5i(jV) + g3 + C2 " 



(-1) 



N 



N(N+l) ' (N+l) 2 ' 

^ + l)(A. + 2)( 2 ^ + 45l - 2W+25 - 3W+25l(Ar) ^ 

2& , 4S_ 2 (A0\ £ 2 -S 2 (N) 3N 2 +N-S 

+ (N + 2) 2 + Af(Af+l) 4 (Af + 2) 2 ' 



(A.4) 



N+\ + ' N + 



f 2 ) 



(A.5) 



^ CI " 25-2 (AO C2 + S 2 (AO C2 - 35 4 (AO - 2S_ 2 (A0 2 - S 2 (N) 2 



|( 1) ^ C(2) + 25- 2 (jV) | 5 2 (iV) 



Af 2 



A/" 2 ' 



(A.6) 
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